Theory of Anisotropic Hopping Transport due to Spiral Correlations in the Spin-Glass 

Phase of Underdoped Cuprates 



in 

O 

o 

^— > 

o 
O 



•4-* 



o 
o 



> 
o 

O 

in 
o 
-i— > 



T3 
O 

o 



X 



Valeri N. Kotov x >0 and Olcg P. Sushkov 2 >Q 

institute of Theoretical Physics, Swiss Federal Institute of Technology (EPFL), 
CH-1015 Lausanne, Switzerland 
2 School of Physics, University of New South Wales, Sydney 2052, Australia 

We study the in-plane resistivity anisotropy in the spin-glass phase of the high-T c cuprates, on 
the basis of holes moving in a spiral spin background. This picture follows from analysis of the 
extended t — J model with Coulomb impurities. In the variable- range hopping regime the resistivity 
anisotropy is found to have a maximum value of around 90%, and it decreases with temperature, in 
excellent agreement with experiments in L^-xSr^CuCU. In our approach the transport anisotropy 
is due to the non-collinearity of the spiral spin state, rather than an intrinsic tendency of the charges 
to self-organize. 



I. INTRODUCTION 

One of the central issues in the physics of the 
high-temperature superconductors is the nature of the 
ground state at low doping. In particular, the possi- 
ble co-existence of ordering tendencies in the spin and 
charge sectors at low temperature is currently being ac- 
tively investigated^* 2 *^ Experimentally at low tempera- 
ture La 2 _ :c Sr :c Cu04 (LSCO) co-doped with Nd (LNSCO) 
exhibits static lattice deformation^ as well as incommen- 
surate (IC) magnetic order at doping x « 1/8. The lat- 
tice deformation indicates the presence of static charge 
order (stripes). Recently similar behavior has also been 
found in La2_ x Ba K Cu04 (LBCO), 5 x sw 1/8. However 
for x < 1/8 and in Nd-free LSCO, charge order has 
not been detected, while IC magnetism persists down to 
x 0.02i&i Thus while IC magnetic order is generically 
observed in the underdoped regime, static charge order 
seems to be confined to the neighborhood of x = 1/8. 
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FIG. 1: (a) (1,1) Spiral spin state, (b) Localized (Coulomb 
trapped) holes, with positions pinned by the random Sr im- 
purities (schematic). The ellipses reflect the symmetry of the 
wave- function ©. The orthorhombic a and b directions are 
shown. 

Theoretically the IC magnetism in the cuprates is of- 
ten modeled as originating from static charge stripes, 
and since these are not universally present, one is forced 
to introduce additional concepts, such as "fluctuating" 
stripes, electronic liquid crystals, etcA^ The microscopic 



origin of such states is still not clear and is currently 
being debated. An alternative explanation for the IC 
magnetism which follows naturally from the t— J model, 
is the formation of a non-collinear spiral state, partially 
relieving the frustration due to the hole motion^ Such 
a state would not co-exist with charge order, as can be 
shown in the context of effective Landau theory^ How- 
ever it is still possible that in the presence of anisotropic 
Dzyaloshinski-Moriya interactions the spin spiral could 
cause weak lattice modulation with period generically 
half that of the spiral. While the spiral ground state 
has a tendency to be unstable toward phase separation, 
we have shown recently that in the extended t—t'—t"—J 
model the spiral can be stabilized by the presence of the 
additional hoppings^i The spiral description was suc- 
cessfully applied to explain magnetic properties of LSCO, 
such as the location of the elastic neutron scattering 
peaks and the change of the incommensurability direc- 
tion by 45° across the superconductor-insulator bound- 
ary [x = 0.055)>i2 The spiral was also proposed as a can- 
didate for the ground state in the spin-glass (insulating) 
phase for 0.02 < x < 0.055I 12 ! 1 ?' 14 

In the present work we address transport properties 
within the spiral framework. Our main motivation comes 
from the recent experimental data in the spin-glass phase 
of LSCO showing transport anisotropies as large as 50% 
to 70%, both in DC and AC measurementsii£*i& It has 
been argued that these data provide indirect support to 
the notion of fluctuating stripes or electronic liquid crys- 
tals, but no quantitative theory exists that takes these 
concepts into account £ Recent infrared experiments that 
carefully identify phonon modes in LSCO once again find 
no charge ordering tendencies in the spin-glass phase of 
this material^ Therefore we take the point of view that 
the ground state has a spiral spin structure, and calculate 
the transport anisotropy in the variable-range hopping 
(VRH) regime for low temperature and frequency, where 
the anisotropy has the largest value. We show that the 
spatial anisotropy of the hole wave-function in a spiral 
translates into anisotropy of the hopping transport. Our 
main result is that for microscopic parameters appropri- 
ate for LSCO (within the t — t' — t" — J model) , the mag- 



2 



nitude of the anisotropy is large (40%-90%, depending on 
temperature) , and close to the one found experimentally. 
Thus we demonstrate that the transport anisotropy data 
can be explained quantitatively within the spiral theory 
which does not involve any tendency of the holes to self- 
organize into charge stripes. The anisotropy was also 
analyzed in Ref. [Tj on the basis of topological defect 
scattering in a spiral; however the results are applicable 
to the quasi-metallic (higher temperature) regime where 
the anisotropy is very small, of the order of several per- 
cent. We emphasize that in the present work we consider 
the low-temperature, strongly-localized regime where it 
is clear from experiment that the anisotropy is the largest 
(10-20 times larger than the value in the metallic regime). 
Finally we mention that the problem was also addressed 
within the spin-charge separation scenario which provides 
a description of the pseudogap phased and could lead to 
effective "insulating" behavior; however within this ap- 
proach the transport anisotropy is linked to the magnetic 
correlation length anisotropy which can be taken from 
experiment but cannot be explicitly calculated. 



A. Summary of Previous Results and Notation 

Our starting point is the description of the spin-glass 
phase (0.02 < x < 0.055) developed in Ref. [l3 which 
is briefly summarized below. First, a single hole resides 
near the points k = (±7r/2, ±7r/2), with a quadratic dis- 
persion around them £k S3 ^-kf + 4ffcf, where k\ is per- 
pendicular to the face of the magnetic Brillouin zone and 
&2 is parallel to the face. Within the effective t—t'—t"—J 
model for LSCO, the parameters are taken to be: tj J — 
3.1, t'/J S3 —0.5, t" / J S3 0.3, where J S3 125 meV. From 
now on we measure energies in units of J (i.e. set J = 1) 
and lengths in units of the lattice spacing a (we set 
a = 1). By using the self-consistent Born approxima- 
tion one finds that these parameters lead to an almost 
isotropic dispersion fi\ S3 /3 2 = ^ 2.2, and a quasipar- 
ticle residue around the nodal points Z S3 0.34. A hole 
trapped by the Coulomb field of the Sr ion generates 
a spiral distortion of the Neel background, which can be 
parametrized as: |i) = e i( >( r ') n W2| |) ; |j) _ e ie(rj)m-<7/2| 
i G "up" sublattice, j G "down" sublattice, where m is 
an arbitrary unit vector perpendicular to the quantiza- 
tion axis of the states | t)> I 4)- The angle #(r), measuring 
deviations from collinearity is given by\H 

s(r) ^^ [l - e " 1 " r(1 + 2CT)1 ' (1) 

where p s ss 0.18 is the spin stiffness, and 1/k ^ 3 — 5 is 
the localization length of the orbital part of the wave 
function. The value of 1/k is extracted from exper- 
imental data^ and is generally expected to increase 
with doping. The unit vector e± — -4^(1, ±1), in the 

usual square-lattice coordinate system. At finite dop- 
ing (0.02 < x < 0.055) the interaction of the long-range 



dipole distortions from holes trapped at different Sr ions 
leads to spiral magnetic order ir^ characterized by average 
9 cx e± • r, as shown in Fig.^a). In a perfect square lat- 
tice the spiral can be directed along any diagonal (1, ±1). 
However from the location of the elastic neutron scatter- 
ing peaksii the incommensurability is determined to be 
along the orthorhombic b direction. In our picture this 
means that the orthorhombic deformation pins the di- 
rection of the spiral, as shown in Fig. ^ and thus we set 
e = e+ = -75 (1) 1)- In what follows the exact nature of 
the pinning mechanism, to be discussed elsewhere, will 
not be important. 

Experimentally the magnetic correlation length is fi- 
nite £ ~ 6 — 20 (decreases with increasing doping) On 
the theoretical side it was argued that £ is finite due to 
topological defects that lead to frustration of the long- 
range spiral orderiiSiH Since the localization length is 
less than the magnetic scale, 1/k < £, we expect that 
the effective description of transport properties, devel- 
oped below in terms of the one-hole wave function, is 
quantitatively valid also at finite (small) doping (i.e. the 
topological frustration mechanism does not affect our 
considerations). It should be noted that the above in- 
equality while being explicit is also the strongest (most 
restrictive) condition we could give. Its refinement would 
have to come out of a detailed theory of IC magnetism 
in the spin-glass phase (i.e. a self-consistent theory that 
takes into account both the effective disorder and the 
spiral formation, generated by the doped holes). Purely 
theoretical arguments aside, transport measurements in 
the doping range 0%-4% suggest that the system is in 
a strongly-localized regime at low temperatures, and we 
thus expect our calculations to be valid as long as dop- 
ing is not too close to the insulator-metal boundary (at 
5.5%). 

The rest of the paper is organized as follows. In Section 
II we analyze the properties of a localized hole in a spiral 
background, and in Section III we use the hole's wave 
function to calculate the in-plane transport anisotropics 
in the variable-range hopping regime. Section IV con- 
tains our conclusions. 



II. LOCALIZED HOLE IN THE PRESENCE OF 
MAGNETIC SPIRAL CORRELATIONS 

The coupling Hsp between the spin of the magnetic 
background (angle 6*(r)) and the orbital wave- function of 
the hole %(r) generates the spiral and has the formtiS 

H SP = -V2Zt J (e ■ X79) X 2 {r)d 2 r . (2) 

The effective Schrodinger equation is then: 

("2 y2 ~ 7 " ^ Zt{e ' w) ) x(r) = ex(r) ■ (3) 

Here the Coulomb potential of the Sr ion which keeps the 
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hole localized is ^- = ^7 , where go is the unit charge and 
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£ e is the effective diefectric constant known to be quite 
farge in the copper oxides, £ e ~ 30 — 100 (increases with 
doping) 4^ The last term in determines the anisotropy 



III. TRANSPORT ANISOTROPIES IN THE 
VARIABLE-RANGE HOPPING (VRH) REGIME 

A. In-plane DC Resistivity Anisotropy 




FIG. 2: Solution of in terms of the expansion @, for A = 
1. Only the first two (largest) components are shown. The 
normalization of xo is arbitrary. Inset: The ratio — Xz/Xo- 



of x( r ) an d is given by 

7 2 t 2 1 
V2Zt(e • V0) J 



\f{nr) - g(K,r) cos(2</?)] , 



j(z) = 1 - (1 + 2z + 2z 2 )e- 2z , f{z) = 2z 2 e 



(4) 



The coordinate system is chosen so that e is parallel to 
the b axis and ip is the polar angle (Fig. The value 
of k (and consequently e) can be found by a variational 
minimization of the total energy^ 
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We will use K ~ 0.3 — 0.4, consistent with experiment^ 
rather than rely on © which contains the uncertainty 
related to the exact value of the dielectric constant. The 
dimensionless parameter A, as defined in JSJ), character- 
izes the coupling between the orbital motion of the hole 
and the deformation of the spin background. Substitut- 
ing the values of Z, (3 appropriate for LSCO (discussed in 
Section LA), we obtain Ass 1, which will be used from 
now on. 

Let us look for solution of © in terms of the series in 
angular harmonics 

x(Kr,tp) = xo(nr) + X2(nr)cos(2(p) + X4:(Kr)cos(4(p) + . . . 

(6) 

Then Eq. Q leads to a set of coupled equations for 
Xi,i — 0,2,4,... We have solved these equations nu- 
merically and the results for the first two harmonics are 
presented in Fig. El The next harmonics are found to be 
small, for example Xi/Xo < 2 x 10 -2 , and can be in fact 
safely neglected. At large distances all the wave- functions 
behave as Xi ^ e _Kr , while at finite distances this behav- 
ior is substantially modified, as shown in Fig. [21 



Below characteristic temperature Tvrh ~ 20 — 30 K, 
the low-temperature resistivity of LSCO in the spin-glass 
phase can be described by the 2D version of the Mott 
VRH formula p ~ exp (T /r) 1/3 i 19 i 20 i 21 i 22 Here T de- 
pends somewhat on doping and sample quality and gen- 
erally decreases when doping increases (and thus con- 
duction becomes easier). As an estimate, for exam- 
ple at 4% doping the data of Ref. |2(J are well fit with 
To ps 500 K; 21 analyzing the curves of Ref. 22 we obtain 
T ps 200 - 300 K, and from Ref. [H we have extracted 
T re 200 K. 

Within the spiral framework the physical idea behind 
the resistivity anisotropy is quite simple. According to 
(|rj|) the hole wave function acquires an elliptic deforma- 
tion induced by the spin spiral, as shown schematically 
in Fig. Bib). Hence hopping in the 6-direction is less 
probable than hopping in the a-direction. Using © , and 
essentially following the derivation of Mott's result^ we 
obtain for the resistivity anisotropy: 



Pb(T) _ Csm 2 (<p)<!> 2 ( K r T ,ip)dip 



Pa{T) 

<I>(Kr, ip) 



L n cos 2 (ip)(fr 2 {nrx , ip)dip 
_ x(kt, ip) 



(7) 



Xo(nr) 

Here r T = ^(To/T) 1 / 3 is the VRH length^ and in or- 
der for the approach to be justified we must certainly 
have Kr-j- > 1. Eq. Q reflects the difference in the wave- 
function overlap for the two electric field directions. The 
(5) factors sm 2 (ip), cos 2 (ip) arise from the fact that the con- 
ductivity varies with the square of the carrier jump dis- 
tance projection in the electric field direction*^ 3 Since 
the hole's dispersion is isotropic ((3\ — P2 — f3), the ex- 
ponential part of the wave function is isotropic and con- 
sequently the resistivity anisotropy is expected to arise 
from the anisotropy of the exponential prefactors, whose 
overlap leads to Q. Both the VRH length and To thus 
remain isotropic within this formulation. 

It is clear from J7J) that as T decreases (i.e. ktt in- 
creases) the anisotropy grows, due to the increase of 
IX2/X0I (Fig. Hpnset)) and consequently the more pro- 
nounced angular dependence of $(kxt,¥>). The results 
are summarized in Fig. |3 where also the evolution of the 
zero temperature anisotropy pb{T = 0)/p a (T = 0) as a 
function of A is shown for completeness in the inset. At 
very low temperature (T < 1 K^i i.e. T/T < 10" 3 ) 
we would have to take into account a crossover to the 
Coulomb gap regime, which however would practically 
not influence the curve in Fig. |3J 

The data of Ref. ^| were taken at temperatures T > 
10 K, meaning that the lowest ratio T/Tq ~ 0.05 (we 
take To ~ 200 K). As the temperature increases be- 
yond Tvrh ~ 20 — 30 K when ktt ~ 1, the approach 
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FIG. 3: In-plane DC resistivity anisotropy in the VRH regime, 
for A = 1. Inset: Maximum (T = 0) anisotropy versus A. 



based on VRH conduction ceases to be valid, as the con- 
duction mechanism changes to impurity band conduction 
and eventually quasi- metallic behavior. In the most rel- 
evant low-temperature range below Tvrh (correspond- 
ing to largest anisotropy), both the calculated magnitude 
of Pb/pa and its temperature dependence are very close 
to the experimental results^ although in these experi- 
ments the temperature is not low enough to penetrate the 
"deep" VRH regime T/Tq < 0.05 where the anisotropy 
should increase even further. It should be also noted that 
there exists quite a bit of uncertainty in the determina- 
tion of To and hence in the determination of the exact 
value of Pb/pa from Fig. [5] 



B. AC Resistivity Anisotropy 



The generalization of the Mott-Shklovskii-Efros for- 
mula to the anisotropic case is straightforward, as it 
amounts to taking into account the non-exponential 
(angle-dependent) part of the wave- function, leading to 
the replacement (with logarithmic accuracy): ln(l/fi) — > 
ln($(Kr^,(p)/r2) = Lv iV , where $ is defined in Q. Tak- 
ing also into account the expressions for e. k we ob- 
tain for the resistivity anisotropy: 
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The expression (jSJ) is plotted in Fig. Due to the log- 
arithmic dependence, in the theoretical limit of zero fre- 
quency the anisotropy would vanish (albeit very slowly) , 
i.e. pb{u> — 0)/p a (uj = 0) = 1. However one must keep 
in mind that the above expressions are not valid at arbi- 
trary low frequencies, as we must have u > T. Generally 
we find that the (maximum) AC anisotropy of 30-40% is 
somewhat smaller than the DC anisotropy. 



P /P 




At finite frequency and temperature the calculation of 
the AC conductivity is a very complicated problem. How- 
ever in the "quantum" limit when the frequency u> T, 
the VRH AC conduction is expected to be dominated by 
resonant absorption by singly occupied pairs (without in- 
volvement of phonons), and is usually relevant in doped 
semiconductors^ The VRH length in the AC regime 
(neglecting for a moment the angular dependence of the 
states) is logarithmic^ r w = (1/k) ln(2|e|/u;), which is 
the main difference from the DC case. From now on we 
denote f2 = w/(2|e|). The formula for r w follows from the 
fact that upon evaluation of the conductivity, the most 
effective pairs are the ones that satisfy^ ui = 21 (r), 
where I(r) = Ioe~ Kr is the overlap integral, and a typ- 
ical estimate of the pref actor is Iq |e|. The func- 
tional dependence of the conductivity in 2D at low fre- 
quency ln(l/51) 3> 1, in the leading logarithmic approxi- 
mation, is given by the Mott-Shklovskii-Efros expression: 
<t{oj) ~ ujr^j[uj + q 2 /r u ], where the second term takes into 
account the Coulomb interaction in the resonant pair^a 



FIG. 4: AC VRH resistivity anisotropy versus Q = w/(2|e|) 
in the regime T < u, and ln(l/fi) » 1, for A = 1. 

When attempting to compare our results to the ex- 
perimental data of Ref. ^| we realize that u> and T are 
not sufficiently low, nor is the difference between them 
sufficiently high to justify the separation into "quan- 
tum" and "thermal" VRH conduction and consequently 
the Mott-Shklovskii-Efros approach. The temperature 
used is T = 13 K and at the lowest frequency of 
uj = 20 cm -1 w 2.5 meV, estimating for the hole energy 
|e| = (3k 2 /2 w 5-10 meV, we have SI ~ 0.13-0.25. More- 
over, upon increasing the frequency to 80 — 100 cm -1 , a 
change of behavior from insulating to conducting takes 
place, a broad peak develops^ and consequently it is 
not clear to us that the AC data are ever in a clean VRH 
regime. In spite of all this the magnitude of the calcu- 
lated anisotropy agrees reasonably well with experiment. 
It should be experimentally possible to lower the tem- 
perature (towards the mK range) as well as uj in order to 
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probe the VRH AC anisotropy more reliably and com- 
pare with our theory. 

IV. DISCUSSION AND CONCLUSIONS 

We would like to reiterate that in our theory there is 
no tendency of the holes to form charge stripe-like struc- 
tures. Nevertheless the resistivity shows anisotropic be- 
havior (and so does the uniform magnetic susceptibil- 
ity, to be discussed separately). It is particularly impor- 
tant, in our view, that the calculations presented in this 
work produce quantitative results, and we would thus 
hope that the successes of the present formulation would 
stimulate further exploration of the spiral and similar 
scenarios, not involving any charge-ordering tendencies. 
The present work is relevant to the spin-glass regime of 
La 2 _ 2; Sr : ;.Cu04, 0.02 < x < 0.055. However we have 
also arguedii*i£ that the spin spiral structure could ex- 
ist in the superconducting state x > 0.055. The the- 
ory naturally explains why the incommensurate direc- 
tion, determined from clastic neutron scattering, rotates 
by 45° exactly at the insulator-superconductor transition 
point. On the other hand it is usually argued that at 
1 = 1/8 both in charge-ordered LNSCOi and LBCO^ 
non-collinear spiral order is not consistent with experi- 
ment. We would indeed not claim that the spiral ground 



state is stable at that particular doping since in fact the 
spiral becomes commensurate with the lattice and there- 
fore significant changes in the ground state could occur 
due to the spin-lattice coupling. 

In conclusion, starting from a spiral ground state which 
unambiguously follows from the extended t — J model, 
we have calculated the in-plane anisotropy of the low- 
temperature DC and low-frequency AC resistivity of 
La2_ 3; Sr 3 ;Cu04 in the spin-glass phase. The theory has 
no fitting parameters and the calculations are performed 
for the variable-range hopping regime. Within the spi- 
ral description the transport anisotropy is due to the 
anisotropy of the hole wave-function since the hole hops 
in a non-collinear (spiral) spin background. The AC 
anisotropy reasonably agrees with experiment in spite of 
the fact that the data are on the border of the VRH 
regime. The experimental data for the DC resistivity 
are well within the VRH regime and here both the cal- 
culated magnitude and temperature dependence of the 
anisotropy agree with experiment very well. 
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